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Abstract. - It has recently been shown that amorphization and melting of ice [|l] were in- 
timately linked. In this letter, we infer from molecular dynamics simulations on the Si02 
system that the extension of the a— quartz melting line in the metastable pressure-temperature 
domain is the pressure-induced amorphization line. It seems therefore likely that melting is the 
physical phenomenon responsible for pressure induced amorphization. Moreover, we show that 
the structure of a "pressure glass" is similar to that of a very rapidly (10^^ to lO" kelvins per 
second) quenched thermal glass. 



Introduction. — 

Pressure induced amorphization (PIA) is an intriguing phenomenon. The first report of this 
transition goes back to 1984 when Mishima et al. |^ found that ice-Ih cooled down to 77 K and 
pressurized to 1 GPa failed to transform into crystalline Ice-IX and became amorphous. Since 
then, a large number of materials were found to undergo similar transitions |^ at different 
pressures and temperatures. In particular, quartz, one of the most studied minerals, was 
found to amorphize at room temperature at pressures between 15 and 25 GPa. A series 
of experiments ||, |^ and molecular dynamics (MD) simulations [0, ||, D have shown that this 
transition is accompanied by an increase in silicon coordination. However, the thermodynamic 
nature of this transition remains very poorly defined. It has been argued on the basis of 
lattice dynamics calculations using the quasi-harmonic approximation and the Born stability 
criteria that the PIA of a— quartz took place as the result of an clastic instability |11 



occurring at approximately 22 GPa. More recently, it has been proposed that a dynamical 
instability occured at lower pressure (18.5 GPa) than the elastic instability, and that the 
sizes and geometries of the simulation boxes dramatically affected the existence and thresholds 
of the reported instabilities. Results obtained by ab initio calculations within the local density 
approximation of the density functional theory showed p3| that the elastic instability occured 
at 29 GPa, but all recent MD calculations |8| Q suggest that amorphization occurs at 
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22 GPa, independently of the size or geometry of the simulation box, so that the role of the 
instability in the pressure induced amorphization process is still a matter of debate | p2[ . 

Hemley et al. ||T^ proposed that PIA of quartz represents the metastable extension of 
the melting curve of Si02- Richet using thermodynamic arguments has inferred that 
the vitrification upon room temperature decompression or upon room pressure heating of 
some high pressure silicate minerals represents crystal fusion. In the first report on the PIA 
of ice, Mishima proposed that the transformation could be related to a metastable melting 
process, and the author recently showed [Q that the pressure-induced transformation towards 
amorphous ice (high density amorphous phase) can be described as a melting transition 
to an extremely unrelaxed amorphous state. This state then slowly relaxes towards an 
energetically more favourable glass-like state. It has been shown using Raman spectroscopy 
that the amorphization induced by heating of stishovite leads to a highly unrelaxed glassy 
state ^ ^1 exhibiting the features of the pressure glass obtained upon room temperature 
compression and decompression of quartz |(^ which then gradually transforms into a relaxed 
glassy state indistinguishable (at least from the vibrational point of view) from the usual 
thermal glass, thus indicating that the two "glassy" states probably belong to the same basin 
in configurational space. 

In this letter, we present results from MD simulations that provide new insights into the 
relationship between PIA and melting, and report the structural and thermodynamic similar- 
ities between amorphous Si02 prepared by static compression of quartz at room temperature 
and quenched silica glass. 

Results and Discussion. - 

MD simulations have been used to investigate the melting and amorphization curves of 
a— quartz in the metastable P-T domain in which the nucleation of the high pressure poly- 
morphs cannot be bypassed in laboratory experiments. Also, the cooling rates used for quench- 
ing liquid Si02 in MD simulations are by far larger than those used by current experimental 
techniques. MD calculations were carried out in an (N,P,T) isothermal-isobaric ensemble, 
for a system of 480 Si02 molecules. The BKS pairwise interaction potential which has 
proven to provide a very accurate description of silica at high pressure was used. A 

system in the a— quartz structure was subjected to compression at different temperatures, or to 
isobaric heating, until either amorphization or melting was observed. In the low-temperature 
high-pressure part of the phase diagram relative to non-equilibrium melting, these amorphous 
states were detected by a combined analysis of the radial distribution function and the energy 
jump at the transition, along with an analysis of the average positions of the particles integrated 
on the duration of the simulation. In the high-temperature low-pressure domain relative to 
the equilibrium melting curve, and thanks to the power of modern parallel computers, we 
could achieve a real thermodynamic determination of the melting points within reasonable 
computation time, which was impossible to carry out only a few years ago on such viscous 
systems . Instead of calculating the differences of free enthalpy between the solid and liquid 
phases, we used an equivalent technique that consisted in performing two phase molecular 
dynamics simulations, and monitored the energy variation of the system with time, along 
with the relative increase or decrease in one of the two phases in the simulation box. With 
computation times of 0.3 ns, it was possible to determine the theoretical melting curve for the 
model, with an accuracy of the order of 50 K at 10 and 5 GPa. 

The corresponding points are reported in the P-T phase diagram of figure 0. They represent 
the numerical coexistence curve between the crystal, and the melt. The term melt is used here 
for both the high temperature liquid state and the high pressure disordered state. Comparing 
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these points with the experimental resuhs for the melting of quartz, it seems clear that 
the amorphization curve is a reasonable extrapolation of the melting curve, given that the 
amorphization process is a continuous extrapolation of the melting process and falls along the 
extrapolated melting line of coesite, the highest pressure form of silica consisting of four-fold 
coordinated silicon clusters. This pressure-induced amorphization line is impossible to measure 
experimentally, because the amorphization process is impaired by the formation of coesite and 
stishovite. It represents the metastable (spinodal) stability limit of the quartz phase. It can 
furthermore be noted that the equilibrium melting curve of quartz in the metastable domain 
{i.e. the coesite stability field) is very close to the experimental coesite melting curve. This 
is due to the fact that quartz and coesite are very similar as far as their local structure 
is concerned (Si04 polyhedra), and even though the solid phases have quite different free 
enthalpies at low temperature, the latter is smoothed at high temperature as the system 
approach the solid-liquid transition curve. Proof of that is the abscence of any break in the 
experimental melting curve of quartz and coesite at the quartz-coesite-liquid triple point. 

Calculation of the metastable extension of the equilibrium melting curve by extrtapolation 
of thermodynamic data JlSf have improved with the measurement of more accurate data at high 
temperature and at high pressure for the Si02 system [ p^ . The recent calculation by Zhang et 
al. [ p8| using the most recent thermodynamic data is quite close to the results obtained from 
the numerical simulation. Although the extrapolated pressures are lower that those in the 
simulation, these authors have argued that these curves should be shifted to higher pressures 
if corrected for the configurational entropy term arising from the distribution of four- five- and 
six-coordinated silicon ions in the melt at high pressure, and for the volume difference brought 
by the appearance of these higher density clusters. In fact, the extrapolation of thermodynamic 
data does not include these higher coordination species, because such a transformation does 
not occur in the transition from quartz to coesite. 

The continuity between PIA and melting should be reflected into some similarities between 
the structural and dynamical properties of the pressure amorphized material (thereafter called 
pressure glass) and of quenched liquids (thermal glasses). The properties of the latter type of 
systems are known to be sensitive to the quenching rate ]l9t . In order to compare pressure 
glasses and thermal glasses, we use as a reference state a sample of quartz, amorphized at 
22 GPa, and then decompressed to room pressure, leading to a pressure glass with a density 
of /5 = 3.07 g • cm~^. Thermal glasses were produced by cooling a high temperature liquid 
(7000 K) at room pressure and constant density p = 2.2 g ■ cm~'^ using three different cooling 
rates 7 = 4.4 ■ 10^^ K • s~\ 7 = 3.5 ■ 10^^ K ■ and 7 = 10^^ K • s^^. These glassy systems 
were then compressed at 300 K to obtain densified glasses with densities of p = 3.07 g • cm~^. 
A second set of thermal glasses was produced by cooling the high temperature liquid (7000 K) 
at a constant density of p = 3.07 g ■ cm"'^, with the same three cooling rates, to avoid the 
final compression step in a highly non-ergodic system. Only the results relative to this last 
set of samples will be discussed below; similar, but less precise results were obtained with the 
samples cooled at p = 2.2 g • cni^^^. 

Static and semi-dynamical quantities such as the coordination number distribution, ring 
statistics (figure ||) and vibrational densities of state (V-DOS, figure ||) were then calculated 
for these systems and compared to the corresponding quantities for the reference Si02 pressure 
glass. 

The coordination distribution shows that the proportion of 4-coordinated Si species de- 
creases with increasing cooling rates, whereas the proportion of 6-fold and especially 5-fold Si 
coordination rises. 

The effect of the cooling rate on coordination distribution is consistent with the pressure 
dependence of the coordination distribution pQ| in liquid Si02; as pressure rises in the 0- 
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15 GPa range, the proportion of 4-fold coordinated Si species decreases whereas that of 5-fold 
(and in a lesser extent to 6-fold) coordinated Si species increases. The higher the cooling rate, 
the closer are the configurations of the supercooled and high temperature liquid (higher " Active 
temperature"). The high temperature liquid at these densities is in a high pressure state, which 
is maintained during the quench resulting in the formation of a strongly frustrated amorphous 
state, because the application of such a high cooling rate "freezes" the system in such a short 
time compared with the mechanical relaxation time that the transformation is quasi-isochoric. 
These two observations represent two different pictures of the same phenomenon. 

On the other hand, the ring statistics show that a high cooling rate reduces the proportion 
of the 5 and 6-membered rings characteristic of a low density structure in favour of small 
chains of 3-membered and 4-membered rings. 

It has been proposed [^l], ^ that the proportion of 3-membered and 4-membered rings 
increases with pressure in densified Si02 glass or pressure-induced amorphous quartz, on the 
basis of the concomitant intensity change of the Raman Di and D2 defect lines associated with 
these rings |2^. Our results for the cooling rate dependence of the ring statistics shows 
that the proportion of these clusters is also enhanced by higher cooling rates, showing that the 
maximum pressure (before decompression) attained during room temperature compression of 
quartz plays a role similar to that of the cooling rate during the quenching of a melt at room 
pressure. 

The coordination distribution in the pressure glass is close to that of a thermal glass obtained 
with a cooling rate 7 ~ 7 • lO^'^ K • s~^ (fig. |^). The ring distribution displays the same 
characteristics, except for 6-membered rings whose proportions (within a few percent) are 
not very significant, due to the greater error bars in the determination of these larger scale 
structures in non-ergodic systems. 

A similar conclusion can be drawn from the vibrational density of states (V-DOS, fig. ||); 
the high-frequency internal stretching modes for the Si04 units are strongly broadened in the 
pressure glass. This broadening is an increasing function of the cooling rates in the thermal 
glass. Another important feature is the peak at 900 cm~^; the glasses with higher cooling rates 
have a narrower peak followed by a plateau, whereas glasses produced by slow cooling of the 
liquid show a broader peak followed by a drop in the V-DOS. Once again, the pressure glass 
exhibits a behaviour intermediate between that observed for the two rapidly cooled systems, 
and is very close to the glass cooled at 7 = 3.5 • 10^'^ K • s~^. These results imply that 
the pressure glass has similar structural and thermodynamic properties as a rapidly cooled 
(10^'^ K • s~^ to 10^^ K • s~^) liquid at the same density and temperature. This supports the 
continuity between the liquid and the pressure glass, and therefore between the melting and 
the amorphization curves. 

A natural question is the relationship between our results and the possibility of a polyamor- 



phic 25, transition from a higher-entropy amorphous form of silica (the pressure glass) 
to a lower-entropy form (the thermal glass). In MD studies, there is no indication of any 
phase transition behaviour during densification, (f^)rp seems negative through the complete 
densification process at room temperature |Q. In order to accurately check the latter propo- 
sition, we carried out very long (N,P,T) simulations (0.134 nanoseconds) on Si02 glass cooled 
at 7 = 4.4 • 10^^ K • s~^ and compressed in the 5-15 GPa range. A histogram of the density 
fluctuations was then plotted and it showed a single yet broad maximum. 

Conclusion. ~ 

All these arguments suggest that below Tg, the metastable extension of the quartz melting 
curve and the amorphization curve are identical, non-equilibrium thermodynamic melting can 
then be considered as the physical process responsible for high-pressure amorphization, at least 
in the case of this theoretical model of silica. 
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Fig. 1. - The quartz- melt coexistence curve (full symbols). This line separates the phase diagram 
into two distinct areas: the quartz metastability field and the region where amorphous silica is more 
stable that quartz. The points represented by circles and squares represent the (vertical) error bars 
on the pressures needed to amorphize the system at constant temperature. In the same manner, the 
diamonds and triangles are the (horizontal) error bars on the temperatures needed to melt the system 
at constant pressure. It can be noticed that pressurization and heating points give the same results 
in the crossover region. The experimental points (crosses) on the experimental equilibrium quartz, 
coesit and stishovite melting curves (thick black line) are reported and the quartz melting curve 
is then extrapolated (thick dashed line, star symbols) to higher pressures and lower temperatures 
using the measured thermodynamic data for quartz (data by Zhang et al. ||18|). The arrow indicates 
that the extrapolated melting curve must be shifted to higher pressures in order to account for the 
configurational entropy and reduces volume changes due to the presence of five and six-fold coordinated 
silicon species ]lq]. 
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Fig. 2. - The structural properties of the glassy systems aX p — 3.07 g ■ crn^'^ and 300 K. The 
coordination (loft) and ring distributions (right) are reported for the three thermal glasses (filled 
symbols) as a function of their cooling rates. The corresponding quantities for the pressure glass 
(open symbols) are close to what would be obtained for a cooling rate 7 ~ 7 • 10^^ K • s~^. 
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Fig. 3. - Vibrational densities of states for the glassy systems at p = 3.07 g • cm~* and 300 K. 
The three thermal glasses produced using different cooling rates are reported. Noticeable differences 
appear due to non-equilibrium averaging; this pseudo-dynamical quantity is very sensitive to large 
scale fluctuations responsible for frequency shifts in the low frequency region. It can bo noticed that 
the rapidly cooled systems and the pressure glass arc closer apart, because high-frequency stretching 
modes of the tetrahedron are strongly broadened in these cases. The observation of another vibrational 
feature at 900 cm~^ leads to the same conclusion. 



